function av_val=wg_integrate(wg_data,threshold)
    [m,ind]=max(max(wg_data')',[],1);
    mask=wg_data>threshold*m;
    mask([1:ind-3 ind+3:end],:)=false;
    av_val=mean(wg_data(mask));
end